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THE RUNGE-KUTTA EQUATIONS BY QUADRATURE METHODS 


SUMMARY 


This report gives a basically new approach to 
the formulation of the classic Runge-Kutta process. 
Dependence on the tedious Taylor expansions is ob- 
viated by a matrix equation which defines the Runge- 
Kutta equations for any order; furthermore, the ele- 
ments of these matrices are obtained quite simply. 
The method of quadratures is used to determine the 
conditions on the parameters which characterize the 
process since these parameters determine the order 
of accuracy of these functions. 

INTRODUCTION 


We briefly indicate the Runge-Kutta process 
which leads to an algebraic system of equations whose 
solutions give the parameters involved. 

Consider the initial value problem 


It is the usual practice in determining the 
parameters in equation (2) to expand these func- 
tions in Taylor series about the initial point (x 0 , yo) 
so that, to a given order, equation (3) is identical 
with the Taylor expansion 

h 2 J 

y(x 0 +h) = y 0 + hf(x 0 , y 0 ) +— r(x 0 ,y 0 ) 

+y, f n (xo,yo) + • • ■ (4) 

which is the true solution of equation (1) . Here the 
superscripts denote total derivatives with respect to 
x, which in terms of the partial derivatives increases 
in complexity as the order increases. On equating 
the respective coefficients of the partial derivative 
in these expansions to a given order in h, we are led 
to the system of equations which determine the param- 
eters in equation (2) . 

When the order of the Runge-Kutta process is 
reasonably high, the expansions in terms of their 
respective partial derivatives become somewhat 
formidable, 1 


y' =f(x,y) y(x 0 )=y 0 . (1) 

We define the sequence of N quantities 


kj = hf(x 0 , y 0 ) 

k 2 = hf(x 0 + cy 2 h, y 0 + /3 2 l k i) 

k 3 = hf(x 0 + a 3 h, y 0 + foikj + (3 32 k 2 ) ( 2) 

k 4 = hf(x 0 + cv 4 h, y 0 + jS^kj + (3^k 3 ) 

k 5 = hf(x 0 + q: gh, y 0 + fafa + /3 52 k 2 + /3 53 k 3 + /3 54 k 4 ) 


N=i 


k N = hf(x 0 + c* N h, y 0 + 2> Nj k.) , 


and the linear form 
N 

y(x 0 +h) « y 0 + Yj k s 
i = 1 


(3) 


With the method we are about to introduce, 
we shall obviate the need for these tedious Taylor 
expansions in order to obtain the system of equations 
which defines the parameters in equation (2) ; 
we will have instead a procedure which is almost 
automatic and extendible to any order. 


THE QUADRATURE METHOD 


Consider the quadrature formula 2 
b N 

/ f(x)dx*£ to f(x ) (5) 

a i = 1 

1 See Kopal's [1] expansion for the fourth-order 
Runge-Kutta formulation; the magnitude of the 
task for higher orders becomes evident. The work 
of Butcher [2-5] and Shanks [6-7] are generali- 
zations based on the approach by Taylor series 
expansion. 


such that the expansion of this in powers of h agrees 
with the exact solution of equation (1) to a prescribed 
number of terms. 


2 We assume a weight factor of unity in equation (5) ; 
obviously the following analysis could be extended to 
include other weight factors. 



which approximates the integral by a linear combi- 
nation of the values of the function for the arbitrary 
zeros x.(i = 1 , 2, . . . , N) . We say that equation 

(5) has the precision m if it is exact for the succes- 
sive powers f(x) = x r (r = 0, i, . . . , m) ( and, 
therefore, also for all polynomials of degree m) . 
This is equivalent to the m+i equations 

b N 

f x r dx = Yj w * x 1 * , r = 0, 1, . . . , m. (6) 

a i=l 1 1 

Suppose now that equation (1) be reformulat- 
ed as 

x 0 +h 

y(x 0 +h) - y(x 0 ) = / f(x,y(x))dx (7) 

*0 

in which form it includes the initial conditions. Let 
us replace the right member of equation (7) by the 
approximating quadrature sum 

N-i 

h I OJ f(x y(x)) (8) 

i= 0 

where we have set 

x i = x 0 + a. 1 h, i = 1, 2, 3, . . . , N-l ( 9) 


Combining equations ( 7) and ( 8) , 

N-l 

y(x 0 +h)-y(x 0 ) « h £ « , f(x., y(x ) ) (12) 

i = 0 

is a quadrature approximation with the degree of 
precision m if equation (10) holds. 

We continue by making another quadrature 
approximation, this time for y(xj in equation (12) . 

For any term f(x., y.) in (12) we can write 
y'(x.) =f(x., y.) =f(x 0 + or. +1 h, y(x 0 + a.^h) ) . 

(13) 

As in equation (7) , we can replace equation 
(13) by the integral form 

x.=x 0 + a. , h 
, . r i u l+l 

y 'V y °“ J f(x, y(x) ) dx (14) 

x o 

and finally, we can replace the integral by the 
approximate quadrature sum 

y< “i’ ' 


y i = y(x i ) =y< x ° +a i+ i h) ’ 
then for a degree of precision m 

h(coi+u;2 + -** aj j s j) =u 0 

h(co 1 x 0 +a;2 x i + *** 00 jq X jsT _ " u l 

. , m m m 

h(w i x 0 + co 2 x i + * * * X N -1 ^ U m 


where the (3 f s correspond to the weight coefficients 
co in equation ( 8) . For a degree of precision m 
equation (10) must be satisfied with the (3 ’s replac- 
ing the co 1 s and u^ replacing u^ in equation 

(11) where 

u / ..r+1 r+1 

x 0+ a 1+ i h r (x 0 + Q;i+1 h) -x 0 

do) V“ i+ i )= J x dx = 771 ’ 

x 0 

r = 0, 1, .. ., m. (15) 


where, as in equation (6) 


x 0 +h 

u = f x r dx = 
r J 

Xo 


. ,.r+l r + 1 

(x Q +h) -x 0 

r + 1 


r = 0, 1, . . . , m . 


Summarizing the above, we can replace y(x.) 
in equation (12) with the requisite precision, if we 1 
set 


(ID 


i-1 


y(x.)*y 0+ h£^ i+1) j+1 f(x., y(x.) 


(16) 


where 

The m+1 linear equations (10) , as we shall 

see, will give us some of the equations in the Runge- h (0 i+l, 2 + ‘" +/3 1+1, i) =U °(“i+l) 

Kutta process. 


2 



xj = (x 0 +a i+1 h > 3 


h (' > 1*1,1 *» ^l+l, 1+1, 1*1 -l) 

(“<♦*) 


(17) 


m \ 

1+1, i X i-l) 


= ui 


^i+1,1^^ 1+1,2^— + /? 

= u m(“i+l) 
i = 1, 2, . . . , N-l 


must be satisfied to obtain the required precision m; 
where, as we shall see, m = p - 1 with p the order 
of the Runge-Kutta process. 

By expanding equation (16) consecutively for 
i = 1, 2, ... and by using the approximations for 
y(x.) in each successive expansion, we will finally 

obtain in place of equation ( 16) 


i-i 


y(x.) * y 0 + Z o * 1+ i. j+ i k j+ i 


(18) 


When this is substituted into equation ( 12) , 
the functions take the form of the Runge-Kutta ex- 
pressions defined by equation (2) and the numerical 
solution equation (3). 


THE RUNGE-KUTTA EQUATIONS 


It will be convenient to put some of the pre- 
vious equations into a more usable form. Thus, for 
equation (10) we can write 


h Jco 2 • * • “to] 


4-i 


= u., j = 0, 1, . . . , m . ( 19) 


Using ( 9) and expanding by the binomial 
theorem, we have 


[ lffl i+i“!+i---“i+i] 


j x 0 3 1 h 


( 20 ) 


Similarly, we can put equation (11) in the 
matrix form 

4 


V h t 1 2 3 " (j+i)] 


3 x 0 H h 


L h 3 J 

so that we can replace equation (19) by 


( 21 ) 


h|w 1 a>2 


10 0 . . . 0. 

1 OL 2 O'! • • • 

1 o?3 a 3 . . . 


L 1 a N a N’ °4 j 


4 

jx 0 3-1 h 


- [‘ I i-o-id 


V 2 

X 0 

3 x 0 H h 


( 22 ) 


3 = 0, 1 m 

If we equate coefficients in equation (22) , 
we get, for any fixed j, some, and finally (for j = m) , 
all of the equations 


jwi <*>»••• w k ] 


100 ... 0 
, 2 m 

L Oi 2 a 2 ■ • • O' 2 

4 2 m 

1 « 3“3 ••• 


j 2 111 

1 a N a N a N 


n ‘ i . — i— I 

L 2 3 (m+1) J 


(23) 


3 



which replace equation (10) and which constitute a 
part of the Runge-Kutta equations. 2 3 


We can treat equation (17) similarly and write 
10 0 . . . 0 . 


h p i+1, 1^ i+1, 2 •" 


1 a 2 a\ . . . a J 2 


1 a . c?. ... ol 

l i i 


x 0 


j *» 


= U j(“l + l) 


i = 1, 2, . 

j = 0, 1. 


(N-l) (24) 

, m 


and as in (21) for u_., we can write for equation ( 15) 


u j ( °i + i ) 


= h 


1 2 13 

a. . — a. 4 - or 


1 v J+ 1 

- )of. 


■i+1 2 “i+1 3 1 + 1 <j+i) H 

j = 0 } 1, . . . , m . 


x J 

x o 

I j 4 1 h I 


(25) 


Substituting this into equation (24) and equat- 
ing coefficients for ail fixed j's, we can summarize 
our result by replacing equation (17) with 


V 


i+1, 2 


(3, 


i+1, i 


(l 2 1 3+0 

= [i a i + l ’■* (J+l) ^ 1+1 J 


a 2 G? 2 • - 

Of 3 Gf 3 • ■ 


a . a 7 

1 l 


i = 2, 3, 

j-i. 2, 


J 

^ J 


a! 

i 


(26) 


(N-l) 

m 


and 


a . , ■ 

l+l 


i+1, 1*^ i+1, 2 + ' ' ' +/? i+l, i. 


(27) 


i = 1, 2 (N-l), 


2 For the solution of equation (23) , with cj . in terms 

of the a r s, see [8] and [9] . 


the latter being the first equation in equations (17) . 
The equations (27) are always assumed in the Runge- 
Kutta process as basic relations between the param- 
eters in equations (2) . 


We might also require that the equations 

P i+l i2 a * +P l+l f 3 a>+ -- + ^i+i, i“i = 2 “l+l 

i = 2, 3, ..., N - 1 
(28) 


hold. These and similar equations obtained from 
equation (26) , however, are not specifically assumed 
here when deriving the Runge-Kutta equations but 
are sometimes invoked in reducing the number of 
Runge-Kutta equations (for p > 4) when solving for the 
parameters. 


We now come to an important step in deriving 
the classical Runge-Kutta equations. We form the 
following weighted expressions using equations (17) 
(but excluding the first equation, or equations (27) ) : 


N-l 


/. co . . u. (a . , ) , 

H . l+l j l+l 

i=i J 

or, using equations ( 26) 


j = 1, 2, . . . , m 


(29) 


3, N 


= Q 


2, N 


where 


2, N 


0 32 

0 

0 


Of 2 Of 5 

m 

• • ■ 0/ 2 

0 42 

0 43 • • 

^4, N-l 


Gf 3 0t\ 

m 

... az 

_ 8 N2 

^N3 

^N, N-l 


2 m 
°N-1 a N-l ’ ‘ ' “N-l 


A 2 

1 


m + 1 



2 012 

••• (m + 1) 

Of 2 



A 2 

1 


m + 1 



2 

*' (m+1) 

Oi 3 

(30) 


A 2 

1 


m + 1 



_2 “n' 

•• (m + 1) 

“n 


(u) 2 W 3 . • • 

“n* 





3, N 


= (oj 3 u) 4 ... w N ) . 


(31) 


4 



If we now set 
c 3^ = P 32 a 2 
c 4^ = P 42 O' 2 + P 43 3 


°N ^ N2 “ 2 + ^ N 3 “ 3 


’• + P N,N-l ff N-l 


1 = 1. 2, .... m-1 
then matrix equation (30) becomes 


(32) 


3,N 


C3 (l) 

cl*) 


( 1 ) 

"N 


C3 


( 2 ) 

( 2 ) 


(m-D 

(ra-i) 


( 2 ) 


(m-1) 


N 


■[HD HD- i (=irl 


(33) 


where we have used equation ( 23) to evaluate the 
matrix product in the right member of equation ( 30) 
except in the last column where it was not possible 
and which we have suppressed in both sides of (33) 
since they yield higher order terms and are not rel- 
evant to our degree of precision. 

Equation (33) gives us another set of Runge- 
Kutta equations. If we return to equations (29) 
and the subsequent analysis, it should be apparent 
that we should be able to find other relevant Runge- 
Kutta equations if we replaced co in (29) by 


oo . . a . . with n = 1 , 

l+l l+l 


2 , 


m-2. 


Therefore, if we write 


(m-1) 
2, N = 


cj 2 

00 2 Of 2 


N 


n 


(m-l)_ 

3,N 


• 

m-2 

• 

m-2 

|_^2 

<*2 ... 


J 

00$ 




CO 3 

« 3 • • • 


a N 


m-2 


m-2 

oo 3 

o' 3 ... 

W N 

% _ 


(34) 


we should be able to replace the row vectors Q ^ ^ 
and Q* * n equation (30) by the (m-1) -rowed 


matrices in equation ( 34) , 
0 (m-l) c = 


This gives 


3, N 

HD 

5 (D 


ft)... 

(?)... 


m-1 


(±) -i (-*-) 

\ m / m \m+i/ 


— (— 
l - 1 \ m+1 j 


(i r) 


(35)J 

where the substitution of C for the matrix in equation 
(33) is evident. Thus, we have extended the number 
of Runge-Kutta equations by equation (35) ; but others, 
as we shall see, are possible. It might be appropri- 

( m-1) 


ate to refer to the elements of Q 


3, N 


and other simi- 


lar elements which we will form as generalized Runge- 
Kutta weights. 

It should be apparent that if we could find all 
the other pertinent sets of generalized Runge-Kutta 
weights, we should be able to find all the relevant 
Runge-Kutta equations. The procedure for doing this 
is reasonably simple. As we obtain new Runge-Kutta 
equations we must examine each for new generalized 
Runge-Kutta weights which in turn may generate new 
equations. The simplest way to recognize these are 

as the coefficients of a ^ , for appropriate integral 

values of n, of the already known Runge-Kutta equa- 
tions. The two specific examples — for orders four 
and five — which will be given should sufficiently 
illustrate the procedure. 

We might finally summarize our results. The 
Runge-Kutta equations for N evaluations in equations 
(2) , with N appropriately chosen, and order p = m-l 
are given by equations (23) , (27) , (35) and all addi- 
tional Runge-Kutta generalized weight coefficients 


that can be found to replace 


n ^ N in equation ( 30) . 


This will result in an appropriate extension of equa- 
tion ( 35) , wherein the matrix which contains aR the 
Runge-Kutta weight coefficients will, for convenience, 
by referred to as Q . 

It might be noted that we also use equation (28). 

JD _ 1 2 


Then, because of equations (32) , c 


i+1 


= — a. , 

2 i+l 


( i = 2, 3, . . . , ( N-i) ) ; and, as may be seen from 
equation (33) , the number of Runge-Kutta equations is 
reduced. 


I. 


in i mi 


(40) 


FOURTH-ORDER RUNGE-KUTTA EQUATIONS 

We will illustrate the preceding work, using 
our quadrature method, for the particular fourth - 
order Runge-Kutta process and show how to derive 
the equations which define the parameters a . and 



We first replace the right member of equa- 
tion (7) by a quadrature sum involving the four 4 
arbitrary zeros (i = 0, 1, 2, 3) , so that corres- 
ponding to equation (12) , we have 

y(x 0 +h) -y(x 0 ) 

* h[o; 1 f(x 0 , y 0 ) +w 2 f(x 1 ,y 1 ) + u> 3 f(x 2 , y 2 ) + w 4 f(x 3 , y 3 )] 

(36) 

where x. and y. are defined by equation ( 9) . We will 
require a degree of precision of three, not four, in 
equations (10) because of the factor h in (36) . 

We will now replace the y (i = 1, 2, 3) in 

(36) to the degree of precision three in accordance 
with the analysis previously given. Thus, the first 
term in equation (36) which will be affected by quad- 
ratures for y(x.) is f(xj, y t ) , which is to be replaced 
by * 1 

f(x 0 +cv 2 h, y 0 + /3 2 j ki) (37) 

as may be seen from equations (2) . By equation 
(14) 

x 0 + a 2 h 

y(xi)-y 0 = L f(x, y (x) dx ) (38) 

x 0 

and after approximating the integral with the quad- 
rature equation ( 16) , and using k t defined in equa- 
tions ( 2) , 

yi - y 0 = (h p 2 i) f ( x o> yo) = P 21 k i (39) 

where, according to equations (17) , the weight co- 
efficient (h 21 ) must satisfy the following four con- 
ditions if the quadrature is to have a degree of pre- 
cision three 

* 

We are anticipating the fact that for the fourth- 
order process only four arbitrary zeros are nec- 
essary. This equality between the order and the 
number of function evaluations, N, only holds for 
orders up to the fourth. For higher orders this 
relation is complex; Butcher [5] deals with this 
problem. 


(h/? 2 i) x o' = u r (« 2 ) , r = 0, 1, 2, 3 

where u r (a 2 ) is given by equation ( 15) . 

In similar fashion we replace y 2 and y 3 in 
equation (36) by making the respective quadratures 
indicated by equation (16) 

Y 2 -yo= (b/3 31 ) f(x 0 , y 0 ) + (h /3 32 ) f (x u y t ) (41) 

y 3 - yo = ( h 04i) M x o> yo) + (k fi 42 ) f( x i> yi) 

+ (h/? 43 ) f (x 2 , y 2 ) (42) 

where (h /3 ..) may again be interpreted as the weight 
coefficients in the quadratures. 

Equation (39) , which defines y x with the 
requisite precision, can be used to replace y t in 
equation ( 41) . Then 

y 2 - yo » /3 31 hf (x 0 , y 0 ) +/3 32 hf (x,, y 0 +/3 21 kj) 

so that, with k t and k 2 defined as in equations ( 2 ) , 

y 2 * yo + 0 31 ki +0 32 k 2 . (43) 

In a similar fashion, by using both equations 
(39) and (43), we can replace equation (42) with 

y3 ~ yo + /?41 k l +P 42 k 2 + /? 43 k 3 (44) 

where the k’s are defined by equations (2) . This 
procedure, as will be observed, replaces y. in 

f(x., y.) with the same degree of precision as the 

degree of precision in the quadrature (36) . 

Again, in order that the required degree of 
precision be at least of order three, the weights in 
equations (41) and (42) must according to (17) satisfy, 
respectively, the equations 

(h foi) xf + (h 0 32 ) xf = u r (ct 3 ) , r = 0, 1, 2, 3 (45) 
and 

(hj 8 «) x, r + (hp a ) xf + (h £ 43 ) xf = u r (a 4 ), 

r = 0, i, 2, 3 (46) 

where (a 3 ) and u ^( a 4 ) are given by equation (15) 
and x. by equations ( 9) . 
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Our problem is to satisfy the conditions (40) , 
(45) , and (46) , so that the repeated quatratures we 
have made represent the Runge-Kutta process to the 
necessary degree of precision. Rather than use the 
general results we have derived, we will illustrate 
our procedure in some detail. 

CONDITIONS FOR A FOURTH-ORDER 
RUNGE-KUTTA PROCESS 

We have seen that in order that the degree of 
precision represented by the bracketed term in equa- 
tion (36) be at least of order three in h, the four con- 
ditions given by equations ( 10 ) must hold. If we re- 
place x. by equation (9) and u. by equation (11) , equa- 
tions ( 10 ) gives 

GOi +<x> 2 + CJ3 +<x>4 = 1 


Let us now consider the expressions formed 

by 

w 2 u. (a 2 ) + <^3 ^(* 23 ) +W 4 U. (a 4 ), i =* 1,2,3 . (53) 

As an example, when i = 2 , this will give, if 
we use equation (15), 

(Xq h) [W2 Q!2 + Ct, 3 a 3 + w 4 ^4] 

(54) 

+ X 0 h 2 [a) 2 Qf 2 +W 3 Gf 3 +C 04 Qf 4 ] 

+ 1/3 h 3 [a? 2*2 2 + *<^ 3 *2 3 + to 4 G! 4 ] 
and by equations [47] this becomes 


|(xgh) + |x 0 h J + |h 3 ) ( 55 ) 


1 

C 0 2 Q '2 +<^3<23 +<^4«4 = 2 
to 2 a 2 + CO3O3 + total = “ 
( jj 2 a \ +00303 + co 404 = ~ 


(47) 


which is free from the weights to . and parameters a . 
and (3... 1 1 

13 


The corresponding expression for equations 
(53) , if we use equations (40) , (45) and (46) , is 


which may have been anticipated by using equation 
(23) with m = 3. 

Four more equations are necessary for the 
Runge-Kutta process which we will now derive, but 
which for convenience we anticipate here. These are: 


x§ h [to 2 a 2 +cl> 3 q ! 3 +CO 4 O 4 } 

+ 2 x 0 h 2 [ ( co 3 jS 3 2 + & 4 0 42) « 2 + ( to 4 0 43 ) a 3 ] (56) 

+ h 3 [ ( CO 3 0 32 + to 4 042) a 2 + ( w 4 043) *2 3 ] 


(*^3032 +CtJ 4042) ( 22 + (<^4043) *23 = “ (48) 

(*^3 0 32 + w 4042)*22 + ( w 4 0 43)*2 3 = ^®) 

(C0 3 £2 3 ) (0 32*2 2 ) +(C0 4 O?4) +(0 42*22 +0 43*23) = 3 ( 50 ) 

( w 4 043)0 32*22 = 24 ’ * 51 ) 

We now return to the conditions ( 40) , ( 45) and 
(46) , which we shall satisfy in the following manner. 
The equations arising from u 0 (a . ) (i= 2, 3, 4) can 

be satisfied identically by using equation (15) , so that 

*2 2 = 0 21 

*2 3 = 0 31 + 0 32 

*24=0 41 + 042 + 043. ( 52 ) 

These are additional relations which hold for the 
fourth-order Runge-Kutta process and may have been 
found directly from equations (27) . 


so that if we equate the coefficients of equations ( 55 ) 
and ( 56) , we obtain three equations for the fourth - 
order Runge-Kutta formulation — equations ( 48) and 
( 49) , and one which is a repetition of an equation in 
equations ( 47) . No new results will be derived if we 
treat similarly (53) when i = 1, 3. 

It will be observed that if we replace co in 

equations (53) by cj . , the resulting expression 

corresponding to equation (55) will again contain terms 
which are free of parameters; in fact, we are led to 
the new equation (50) . This is an extension of the 
process we have described and which leads to the con- 
cept of the generalized Runge-Kutta weight coefficients; 
and, as we have pointed out, if we can find all sets of 
these coefficients which are relevant to a particular 
order of the Runge-Kutta process, we will produce 
all the Runge- Kutta equations . We continue with our 
coefficients by seeking contributions to these from 
each new equation as it is generated. These are, in 

our present case, the coefficients of a . and a } in 

ii 

equations ( 48) and ( 49) which contribute a new set of 
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generalized Runge-Kutta weight coefficients which 
generate the new equation (5i) 5 . 

We can summarize these results by the ex- 
tension of the matrix equation (35) which was derived 
by a generalization of the elementary concepts pre- 
sented here for the fourth-order case. This gives 


where the c ^ (i = 3, 4, 5, 6) are defined by equa- 
tions (32) and where the new quantities introduced 
are given by 

^32 0 ;• 6 ] 

P 42 P 43 0 


E 


n n 

oj 3 o ( 3 0)40(4 . 




n 

N 


o> 3 cu 4 


0 32 a 2 £ 32^2 

0) 3 Q?3 0)40(4 


(0 42 &2 +P 43 01 3) ((3 42 G? 2 +0 43 ff 1) 

W4/3 43 0 



^N2 ^N3 


■P 


N, N-l J 


(59) 



which represents the results given by equations (48) 
through ( 5i ) . 

THE FIFTH-ORDER RUNGE-KUTTA EQUATIONS 

Several aspects on formulating the matrix £2, 
which contains all the generalized Runge-Kutta weight 
coefficients, can best be illustrated by considering 
the fifth-order case, which required at least six eval- 
uations for a solution [4, 6], 

We start with the applicable Runge-Kutta equa- 
tions given by equations (23) and (27) with N = 6 and 
m = 4. We then enumerate the generalized Runge- 
Kutta weight coefficients; starting initially with the 
three rows given by equations (34) and followed by 
the subsequent extensions as new Runge-Kutta equa- 
tions are derived. This gives 



“ cu 3 

CO 4 

^5 

O >g 


co 3 a 3 

0 ) 40(4 

0 ) 50(5 

o) 6 Q( 6 


Cx) 3 al 

0 ) 40(4 

o >5 Of 5 

0 )$Ol\ 


(o) 

( 0 ) 

( 0 ) 



73 

74 

75 

0 

a = 

(0) 

. <°> 

^ ( 0 ) 



73 «3 

74 Of 4 

75 <*5 

0 


(1) 

(1) 

(1) 



73 

74 

75 

0 


(1) 

(l) 

(1) 

(1) 


CQ 3 C 3 

CO 4 C 4 

^5 C 5 

0)0 Cg 


(74 l 3 43 + 75^°l 3 53> 

75 ( °^54 

0 

0 


— 





5 Other generalized Runge-Kutta weight coefficients 

which may be obtained from equations (48) through 
(51) may be ignored, if we keep in mind that in the 
Runge-Kutta equations we seek no term which is a 
product of the ck’s and (3 's may have the sum of its 

powers exceed m = p - 1. 


r (n) (n) (nfl 

- 72 73 7 n J » n = 0, 1, 

The Runge-Kutta equations for the fifth order 
will then be given by extending equation (35) : 


i (i) i (i) i (i) 



§ [i ft)]* 
I ft (I)]* 
I [I (I)] 


|j [I ft I)] — — _ 

(60) 

where C is a (4x3) matrix defined in equation (33) , 
and where the numerical values of the elements in the 
right hand matrix in equation (60) may be readily 
verified partly by equation (23) and partly by equation 
(30) ; 6 and the other by using the new Runge-Kutta 
equations for the fifth order as they are generated. 

The dashes in this matrix indicate that there are no 
results relevant to our problem. 7 

If two of these equations (those obtained from 
the fifth and sixth row of & in equation (58) and mark- 
ed with asterisks in equation (60) ) are added, instead 

6 The first three rows may be directly obtained from 
equation (35) . 

7 Except, perhaps, for higher order terms which may 
be of interest in evaluating error terms. 


8 



of being satisfied simultaneously, we get the set of 
Runge-Kutta equations for the fifth order customarily 
given in the literature. 8 

THE GENERALIZED RUNGE-KUTTA 
WEIGHT MATRIX 

We shall now indicate a device whereby we 
can generate the matrix of equation [58] quite easily. 
The convenience of this device would be realized if 
we attempted to generate ft from the basic definition 
of the generalized Runge-Kutta weights previously 
given. Even for the relatively simple case of the 
fifth order the latter process is extremely tedious; 
furthermore, the evaluation of the numerical matrix 
[60] is laborious and there is, additionally, the pos- 
sibility of overlooking some rows of ft if the order is 
high. 


need these in this case) ; and let us replace cj ^ as 


we did before, by , y. ^ 


, (o) ( (o)\ (o)/ (i)\ 

by V Vi )' y i ( y i )> 

(Vi ■ 

/ < (° 

[ y > ) = 74 


and designate these 
. , so that for ex- 


ample 

(o) 


(°) (o) (0)„ 

74 P 43 + 75 P 53 + 76 0 63 


) ( 0 ) ( 0 ) 

P 43 “4 + 75 P 53 <*5 + 76 P 63<*6 


where y 6 = 0 in the present case. These again, as 


we will show, are generalized Runge-Kutta weight 
coefficients. In this case, after rejecting those ele- 
ments whose degrees are too high, we retain only 



the elements of the eighth row in (58) . 


The procedure we will now introduce is quite 
simple and we will illustrate it for equation (58) . 

We start with the first three rows as required by 
equations (34) with m = 4. We next replace oj . 

(i = 3, 4, 5, 6 ) consecutively by its corresponding 

terms y.^, y.^ , ... (i=3, 4, 5, 6 ) where 

the latter terms are defined by (59) . This contri- 
butes the fourth and sixth rows to the ft -matrix as 

we reject y ^ , 7 ^^ , ... since the total degree of 

these terms (the sum of the exponents of the a's and 
j 3 ’s) would be too high. 9 We now repeat this procedure 
replacing a>. , as before, in the d. a. and c 0 . a 2 rows 

of equation ( 58) — retaining only those terms within 
the limitations imposed on the total degree of the re- 
sulting elements. Thus, we reject y ^ a 2 whose 

total degree is three and retain only the elements 
( 0 ) 

y. a. in the fifth row. 

1 1 

We continue to make analogous substitutions in 
the elements of the new rows generated. Let us 
(0) (0) 

write y. (to .) for y. to show, as we may see by 
equation (59) , that y|°^ explicity involves to (the 
same is true for y/^, . . . though we will not 

8 The Runge-Kutta equations obtained here may be 
checked with the equations given by Luther and 
Konen [10] who use a slightly different notation. 

9 The total degree ( see footnote 5) of the fifth-order 
generalized Runge-Kutta weight coefficients in 
equations (58) cannot exceed m - 2. This is be- 
cause its products with c.^ 1 ) (whose total degree 

is two) will give the maximum total degree permis- 
sible for the fifth-order case. 


One further addition to our catalogue of terms 

must be made. We should, generally, add to our 

. . , ( 1 ) n (2) n 

matrix such rows as oj.c. a. , w,c. a. ... 

11 1 ii 1 

(i= 3, 4, ... N; n = 0, 1, ... , (m - 4) ) . In our 

specific case only co.c.^ (the seventh row in equation 

(58) ) is tenable. These elements, too, maybe 
subjected to the substitutions indicated above to give 

us terms like y.^c.^, y.^c.^ a., etc. 

r i 1 '1 1 1 


We will now justify the rules we have just 
given for constructing the ft -matrix and show that the 
entities we obtain are generalized Runge-Kutta weights. 
Since the generalization to any order would be obvious 
in what follows, we will, for the sake of brevity, con- 
tinue with the demonstration for the fifth -order case 
with six evaluations. 


Using equation (30) , we can write 


C 03 

0 ) 30^3 

W3«3 


o; 6 

OjQOtQ 


U>2 Ci>3 


Ct)2^2 aj 3 Qf 3 


P 32 0 0 0 

P 42 P 43 0 0 

P 52 P 53 P 54 0 
P 62 P G3 064/3 65 


U> 6 

OJ 6 <*61 


W6<*6 


*<*6 


a 2 <*2 <*2 
<*3 <*§ o>l 

2 3 

a 4 a% <24 

2 2 

<*5 a% as 


1^ 

2 

<*2 

1 3 

3 

1 4 
4 “ 2 

1 

2 

l 3 

1 4 

2 

<*3 

3 013 

4 “ 


“6 


1 2 1 3 1 4 


"«e 


( 61 ) 


and using equation (59) on the left side, we get 
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„ <°> 

(0) 

(0) 

(0) 

72 

73 

74 

75 

(1) 

(1) 

(1) 

(1) 

72 

73 

74 

75 

,, < 2 > 

(2) 

(2) 

(2) 

72 

73 

74 

75 


5 (I) 

■ I (;) 

I (I) 


i (;) 

s (I) 


a 2 a 2 &2 


a 5 oi\ al 


where the right member of equation (61) has been 
evaluated by equation (23). From equation (62) it 

may be seen that since the terms y ^ ot} (n = 0, 1,2; 

j = 0, 1,2,3) are coefficients of in Runge-Kutta 
equations, they fulfill the basic definition of general- 
ized Runge-Kutta weight coefficients. 10 


If now in equation (61) we replace a;, by y. 


we get 1 

r! (°) 


m- 

h- 


( 2 ) / ( 0 ) 


< 0 ) ( ( 0 )\ 

v. h ) 

* (v, W ) 

” (V) 


a 2 q* 2 a 2 


a 3 oil 


o ?4 <24 0?4 


(o) (0) 

72 73 

(°) (0) 

72 Qf 2 73 <*3 

(°) 2 (°) 2 

I 72 73 <^3 


(°) 

76 ^6 

( 0 ) 2 
76 a?g 




' 2 — — ^,4 1 


2 (it) 
i(H) 
j ft t) 


i (a) i (i A' 

3 \4 5/ 4 \5 6/ 

i ( ± i\ — 

3 \5 6/ 


where the last matrix product was evaluated by using 
equation (62). Thus, we have demonstrated that en- 


tities of the type y 




(where m may be 


0) are also generalized Runge-Kutta weight coeffi- 
cients. These simpler cases may be compounded 

10 (n) 
Note from equations (59) that the missing y 6 

( n = 0, 1 , 2) in ( 62) are zero. 

11 (o) r <0,1 (o) r (oil n 

75 7, = 76 Y, =« 


indefinitely to produce coefficients for any other 
order Runge-Kutta process. 


If in taking the product in the left member of 
equation (61) we use equations (32), we get, instead 
of equation (62) , the matrix equation 



which is equation (35) with its right member unex- 
panded. If we expand the left member of equation 

(64) , we see that terms like oj. c/ n ^ a) are also 

l i i 

generalized Runge-Kutta weight coefficients. Evi- 
dently then, we can replace u>. in y. ^ co. by terms 
of this type also . 12 

( 0 ) 

If in equation (64) we replace w. by y . in 

the terms and 6 defined by equations (34) , 
the right member can be evaluated by using equations 
(63) . Thus, in the Runge-Kutta weights of the type 

(n) q (0) 

o^c^ a . we may replace the to. by y^ ; and by an 

extension of the preceding arguments gj. may also be 

replaced by y. ••• , etc., producing terms 

of the type c.^ 01 ^ 

i i i 


It should be recalled that equation (64) (or its 
general form for higher orders) is our basic matrix 
equation, which still gives valid Runge-Kutta equa- 
tions with any generalized weights which may replace 
the elementary weights in and ^2 3, 6 . We should 
keep this in mind below as we evaluate the matrix 
product ftC for the sixth order. 


THE SIXTH-ORDER RUNGE-KUTTA EQUATIONS 

In Table I we have enumerated the generalized 
Runge-Kutta weight coefficients for the sixth-order 
process. We show only the first element of each row 
since the other elements in the row are simply obtained 
by advancing the subscripts to N. The ^-matrix can, 
of course, be easily expanded from the first four rows 
in Table I by the simple rules we have just given. It 
will be observed that if we limit ourselves to weight 
coefficients whose total degree is less than three we 

12 See, for instance, y 3 ^ ^(.0.0. in Table I. 
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I 


will obtain the ft-matrix in equation (58) for the fifth 
order. 13 It is evident, then, that for any higher order 
Runge-Kutta process we need merely extend the list 
in Table I. 

In the last four columns we give the numerical 
values corresponding to the matrix product ft C where 
these are pertinent; 14 and we shall next consider how 
these are evaluated. 


THE EVALUAT ION OF THE NUMER I CAL 
MATRIX 


The numerical values for the first four rows 
in Table I are given by equation (35) ; and we will 
illustrate how, by using previously derived results, 
we can determine the other numerical values shown 
in this table. For the sake of brevity we will limit 
ourselves only to those results just requisite to the 
sixth-order case with N evaluations. The extension 
to higher orders is amply evident. 


From what has already been established, we 
can write the matrix equation, valid for n = 0, 1, 2: 



include those of total degree m - 2. 

14 C is a (N x 4) matrix. Butcher [4] gives a con- 
struction for the sixth order with N = 7; Hilta [3] 
originally gave one with N = 8. 


where we have evaluated the right side by using 
equation (62) or (63). These, and subsequent re- 
sults, have been used in the appropriate place in 
Table I. 


Again, starting with the basic equation 




we again get the numerical results by using equation 
(63). 


Similarly, we can write 
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TABLE I. THE GENERALIZED RUNGE-KUTTA WEIGHT COEFFICIENTS FOR THE SIXTH ORDER 



Total 

Degree 

0 

1 

2 
3 
1 
2 
3 
2 
3 
3 
2 
3 
3 
3 
3 
3 
2 
3 
3 
3 


1 1 
2 3 

1 1 

2 4 

1 l 

2 5 

1 i 

2 6 

1 i i 

2 3 4 

1 1 1 
2 4 5 

1 i i 

2 5 6 

1 i i 

2 3 5 

1 i i 

2 4 6 

1 i i 

2 3 6 

k k k k 

2 3 4 5 

ill! 
2 4 5 6 

1 1 1 l 

2 3 4 6 

1 i i i 

2 3 5 6 

1 1 1 1 

2 3 4 5 

1 i i i 

2 3 2 6 

1 i i 

2 2 5 

ill 
2 2 6 

1 i i* * 

2 3 6 

1 i i i 

2 2 5 6 


1 1 _ 
3 4 

3 5 
1 1 
3 6 


i i i 

3 4 5 
1 1 1 
3 5 6 


i t i 
3 4 6 


1 1 1 1 
3 4 5 6 


1 _ - 
6 


1 1 i * * 

3 2 6 


k k 

4 5 
1 i 
4 6 


i 1 1 

4 5 6 


1 1 
5 6 


* Only the first weight coefficient (element) of each row of the ft -matrix is shown; 
the other elements in the row are simply obtained by advancing the subscripts to N. 
In the last four columns are given the values of ftC. 

** These two are easily seen to be repetitions of each other. 
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( 69 ) 



73 


( 1 ) 


1 

2 





( 68 ) 


3 

N 


where the last evaluation was. made by using equation 
(62) ; so that, combining the results in equations 
(67) and <6 8), we get 
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